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Abstract 

We perform a new, recursive reduction of one-loop n-point rank R tensor Feynman integrals [in short: (n, R) -integrals] for n < 6 
with R < n by representing (n, i?)-integrals in terms of (n, R — 1)- and (n — 1, R — l)-integrals. We use the known representation 
of tensor integrals in terms of scalar integrals in higher dimension, which are then reduced by recurrence relations to integrals in 
generic dimension. With a systematic application of metric tensor representations in terms of chords, and by decomposing and 
recombining these representations, a recursive reduction for the tensors is found. The procedure represents a compact, sequential 
algorithm for numerical evaluations of tensor Feynman integrals appearing in next-to-leading order contributions to massless and 
massive three- and four- particle production at LHC and ILC, as well as at meson factories. 
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1. Introduction 

For the evaluation of next-to-leading order contributions to 
processes at high energy colliders like LHC and ILC, but also 
at meson factories, one needs an efficient and reliable treatment 
of (n, i?)-integrals, i.e. ri-point Feynman integrals with tensor 
rank R. Typically n < 6 and R < n may be needed for final 
states with massive particles. For n < 4 the Passarino-Veltman 
reduction [1] may be applied. For n = 5, 6 there are a variety 
of reduction schemes; for an overview, see e.g. [2, 3, 4]. In 
a recent article [5], we derived such a tensor reduction scheme 
for pentagons and hexagons using the Davydychev-Tarasov ap- 
proach [6, 7] for tensors of rank R < 3. Recurrence relations to 
reduce dimensions and indices have been applied with system- 
atic use of signed minors [8] as described in [9]. Simplifications 
were derived in [5] for n — 5, R — 2, 3 and n = 6, R — 2, 3, 4. 
A corresponding numerical Mathematica package hexagon. m is 
publicly available at http://prac.us.edu.pl/^gluza/hexagon/ [10]. 

In this article, we present a new, recursive tensor reduction 
which extends the reduction up to (6,6)- and (5, 5)-integrals. 
In this reduction tensor integrals with n < 4 occur which can 
also be dealt with in a recursive manner. This will be sufficient 
to evaluate the one-loop amplitudes of four-particle production 
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at LHC and ILC. The new reductions rest on a new master for- 
mula, equation (16) for five-point functions and corresponding 
ones for simpler functions. In principle, the tensor integrals 
with n < 4 might be treated following [1], e.g. with the Fortran 
package LoopTools/FF [11, 12]. LoopTools treats loops with 
massive propagators for n < 5, 1 and Golem [13] with massless 
propagators for n < 6. Unfortunately, there is no publicly avail- 
able numerical Fortran package with a stable treatment of both 
massive and massless particles in the loop. In this situation, it 
appears natural to work out the complete reduction scheme for 
the whole chain of tensors in a systematic way. 

The tensor reductions given for n < 4 in [1] express the 
tensors in terms of scalar 1- to 4-point functions. For tensor 
5-point functions, reductions to tensor 4-point functions with 
rank less by one have been presented in [14] and [15, 16], and 
also in [17, 18] tensor recursions are discussed. A representa- 
tion of scalar (N + l)-point functions (including integrals with 
powers of the loop momentum in the numerator) in terms of TV- 
point functions in TV integer dimensions was derived in [19]. 
For TV = 4, such a representation was derived already in [8]. 
The general case of tensor integrals using dimensional regular- 
ization was treated in a series of papers [20, 21, 22, 23], thereby 
allowing also for massless particles. 

Our reductions express (n, i?)-integrals recursively in terms 
of (n, R — 1)- and (n — 1, R — l)-integrals for n — 2 ... 6. 
Although all approaches have identical basis elements and thus 
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have to have equivalent tensor coefficients when compared after 
complete reduction, we would like to stress that they allow for 
quite different algorithmic realizations. 

The article is organized as follows. Basic formulae are in- 
troduced in Section 2. Section 3 contains the main result, the 
recursive tensor reduction, based on master representations for 
the (5, R) -integrals. As a demonstration, we derive the (5, 4)- 
integrals in more detail in Section 4. Section 5 is a short com- 
ment on the (6, i?)-integral reductions. In Section 6 we state 
some properties of the auxiliary vectors used in the recursions, 
which also allow for an alternative and simple derivation of our 
master formula. We finish with a short Summary. 

2. Basic formulae 

We study Feynman tensor integrals in the generic dimension 
d = 4 — 2e with n external legs: 
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where the denominators Cj have indices Vj and chords qf 

cj = (k — qj) 2 — nij + is. (2) 

We will assume Uj = 1 in the following, but a generalization of 
the results to arbitrary indices is straightforward. 

The iteration of reduction steps will be performed until the 
level of (n, 0)-integrals with n < 4 is reached. In this chain, 
the following well-known scalar reductions [9] are needed: 
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The simplest, but typical tensor is the vector integral: 
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The integrals I„f are scalar n-point integrals, obtained from 
I n by raising the index of line i by one unit (i/j = 2 then) 
and replacing the generic dimension by dimension d + 2. 2 In 
a next step, we apply the recursion relations derived in [7, 9] 
in order to eliminate the shifts of dimension and indices. The 
details of the derivations, which are relatively easy for lower 
rank tensors, get complicated, due to many cancellations, for 
tensors of higher ranks. 

The recursion relation for the vector coefficients in (4) reads 
for n < 5: 
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For n = 1 the second term in (5) vanishes and for n > 5 the 
denominator in (5) is () = 0. The case n = 6 is of practical 



importance and will be discussed below. 7„ is the scalar inte- 
gral, and the integral where line s has been scratched. 3 
The objects like (fj are signed minors, and Q n is the mod- 
ified Cayley determinant. For explicit definitions, see [8] or 
Appendix A of [5]. 4 Thus, we can write for the vector n-point 
function [i.e. the (n, l)-integral] (4) for n < 5 : 
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In (6) we introduced the auxiliary vectors Q^: 
( S ) 



i=l K,r 



-, s = 0, . . . , n. 
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Vectors (7) are universal and will appear in more involved re- 
ductions again. Indeed, (6) is what we want to obtain further on, 
i.e. we will look for analogous relations for higher rank tensors 
in the following. 

Equation (6) is essentially due to recursion relation (30) of 
[9], which reduces simultaneously dimension and index (let us 
call it type I recursion). For (n, i?)-integrals with rank R > 2, 
a complication arises due to the appearence of the -tensor 
in the reduction to scalar functions, as may be seen from the 
simplest case of an (n, 2)-integral: 
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with Vij = 1 + Sij. In d — 4 dimensions, one may eliminate 
g^v by expressing it in terms of the n different chords of the 
integral: 
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For n = 6,5, see [9]. For n < 5, we have to introduce extra 
terms [24], defined with the aid of: 

V* = e^ pa {q 1 -q i ) x {q 2 -q i ) p {q z -q i ) ai (13) 



(14) 



where 



V^V\ = q\q»q v 2 + q 2 2 q^ - ( qi q 2 ) {q^ + 

- blql - (gig 2 ) 2 ] g^- (15) 



2 Analogously, has indices Ui = 2 and Vj 

d + 4 dimensions, etc. 



3 Analogously, in ij,^^ li nes * and t have been scratched, etc. 
2 and is defined in 4 The Gram determinant is \2qjq^\, j, k = 1, — , n — 1, and for q n = 

we have () n = -\2qjq k \. 
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In (15), the q±, qi are short for the two 4-vectors in (14), (gj — 
q$)\i = 1,2. We just mention that in our conventions det (g^) = 
— 1 and e 0123 — — e i23 = +1. Further, v 2 = |()4- It is also 
interesting to note that the contractions of the sums appearing 
in (11) and (12) with their corresponding extra terms vanish, 
i.e. these terms are orthogonal. 

Now, we are ready to derive a systematic recursion algo- 
rithm. 

3. The (5, .Reintegrate 

The scalar (5, 0)-integral is given in (3), and the vector (5, 1)- 
integrals in (6). Applying recursion relations, one may derive 
the following master formula for the (5, i?)-integrals: 
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The formula is a generalization of (6). Equation (16) is given 
implicitly in [5], for n = 5, R < 3: by equation (3.7) for R = 2 
and by equation (3.19) for R = 3, both in combination with 
equations (2.1) and (2.2) (these both for n — 4). For R = 4 see 
also section 4 and, alternatively, see also section 6. 

As a consequence, (16) has a completely new tensor struc- 
ture compared to what one is used to. The explicit evaluation 
of (16) will be discussed below for R = 3,4,5. The (5,2)- 
integrals follow immediately from (16) since the vector inte- 
grals It? and l£' a on the r.h.s. are known from (6). 

3.1. The (5, 3)-integral recursion family 

The master formula (16) for (5, 3)-integrals has on the r.h.s 
(5, 2)- and (4, 2)-integrals: 



(17) 



After eliminating the g^ v in (8) with the aid of (1 1), we obtain: 



-I 



[d+],s 



(20) 



The notations Q s '" and v 3 '" mean that from the five chords of 
the 5-point function the chord s is excluded, such that these vec- 
tors are constructed from four chords - as given in (1 1). In fact, 
(20) together with (19) is the reduction of the (4, 2)-integrals, 
when it is combined with (6) and: 
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where the upper indices s, t in the Q^'-vectors are again intro- 
duced for the scratched lines. Observe that for i = s,t and 
u = s,t there are no contributions so that indeed the indices 
are running only over three values and objects like (°^) 5 can 
indeed be read as (9) , In this way (21) is consistent with (6). 
If one is interested in 4-point functions from the very beginning 
one avoids of course this clumsy notation - but for the present 
purpose of reducing 5-point tensors to scalars it appears ade- 
quate to demonstrate this at least once. 

This completes the (5, 3)-, (4, 2)-, and (3, l)-integral recur- 
sions. 

3.2. The (5, 4) -integral recursion family 

For the higher tensors of the 5-point function we need cor- 
respondingly higher tensors of the 4-point functions, and the 
corresponding extra terms related to the elimination of the g^ v 
have to be derived. 

Thus, as a next step we seek a representation for the (4, 3)- 
integrals which is needed for the (5, 4)-integral recursion: 



The (5, 2)-integrals have already been discussed and may be 
expressed by (5, 1)- and (4, l)-integrals according to (16) and 
(6). For the (4, 2)-integrals, we need now a reduction analogous 
to the master formula (16) and may use as a starting point (8). 
The dimensional shifts in (8) may be treated with a reduction of 
type I: 



[d+] 2 



(0s\ 

\JSJ 



5 jld+hs 



% 3 



4,i 



( tS ) 

\jsJ_5 T [d+],st 
/ j (s\ J 3,i 



5 



W5 T [d+],s 



(18) 



and another one, needed for the reduction of the dimension only 
(not of an index), let us call it reduction of type II. Equation (31) 
of [9] yields: 
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but is also of interest in its own. In the following, for the ease 
of notation, we will drop the scratches of line s. 

A systematic application of recursions of type I results in: 
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In 1% ^ d+ \ besides (known from (19)), also enters. 
It may be reduced by a recursion of type II quite similar to (19): 
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Finally, our representation (23) of 1% vX contains the integrals 
1$ and we have to reduce them also. Application of recursion 
relations yields the analogue of (20) for n = 4: 

i 3 — J 3 2 ^« 3 ?t\ ' ' 

u=l U/4 

We made use of the definition (7), which becomes here: 
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and of the representation (12) in order to express: 
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This completes the (5, 4)-, (4, 3)-, and (3, 2)-integral recur- 
sions. 

3.3. The (5, 5)-integral recursion family 

For the tensor of rank 5 of the 5-point function we need 
further the tensor of rank 4 of the 4-point function. 
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Again, the systematic application of the recursion relations re- 
sults in 



(31) 



with 
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where G^ v and I^ d+ ^ are given in (24) and (25), respectively, 
and 
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The J^" 1 "]'* is known from (26), the from (19), and for 
completeness we specify also the other recursions of type IT. 
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In (31) enters also the j^ l,x < t 
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which we evaluate to be: 

4 



jfiuX,t jpu.t Qt,\ X * jp,v y tu Qt.} 



u=l 



_ Qt,fi\j-v,[d+],t _ Qt,uXjfJ-,[d+],t 

where G*'^ is given in (29) and /^ [d+1 <* in (33). 
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It remains to evaluate the tensor I^' for which we get: 
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The auxiliary vector Q*"'", defined analogously to (28), van- 
ishes when !) = f or f = ti. The sum over w in (37) therefore 
consists of two terms, v = i,i , where i < i and i, i ^ t, u. 

with a 4 - 

basis in (37): 
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Finally, in (37) appear the integrals: 
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and 
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The last equality in (42) follows from: 
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This completes the (5, 5)-, (4, 4)-, (3, 3)-, (2, 2)- and (1, 1)- 
integral recursions. 



4. Derivation of the master formula for (5, 4)-integrals 

As mentioned before, (16) is implicitly contained in [5], but 
only up to R = 3. We will make use of it, however, up to R = 5 
at least. To demonstrate the general approach of how to obtain 
this relation, we give details of its analytic proof for R = 4. In 
particular we observe huge cancellations of higher dimensional 
integrals in (2.4) of [5]. To start with, as usual we write this 
relation as: 
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where (10) has to be used. The highest dimensional integral 

A d +f -,u 

occuring now is nijkil^ ijM , with riijki = VijV ijk Vij kl ,Vij kl = 
1 + 5u + Sji + Ski etc., for which we need the recursion relation: 
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In the sum of the last three terms of (45), to be abbreviated 
as [ijk]( l \ it is understood that there occur no equal indices 
among i,j,k (otherwise it would be contained as a 4-point func- 
tion in the second term on the right-hand side of (45)). With 

the remaining factor VijVijk of the integral nijkll^tjkl we can 
rewrite: 



Us 



., .1)5 Ad+f 
1 '* Q 1 5,ik 



,. (1)5 T jd+) 



(46) 



and as a result we have: 



Vkl Th h > kl 



Ad+f 
() 5 ^ k 



5 

E ■ 

s=l,s^i.j ,k 



1 

Os 



i/ r, 



(i) 



5 r[d+] 3 

-t r 



Or 



(Or 



5,j'i 



- 5 r[«i+] J 



(47) 



where (46) has already completely cancelled against the last 
three terms of the second sum of (2.4) in [5]. With the further 
recursion: 
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one observes that the term cancels against the last row of 

(47) and the remaining 1$ integrals combine to I^^jk according 
to (3.18) of [5] with the result: 
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This proves our statement: The first term corresponds to the 
(5,3)-integral as claimed in (16). The rest, according to (2.3) 
of [5], corresponds to the (4,3)-integral, again corresponding to 
(16). A similar proof for n = 5,R = 5, equation (30), is a 
bit more lengthy, but see also the derivation in section 6. We 
would like to mention that a pedagogical introduction to the 
techniques applied may be found in [25]. 



(45) 5. The (6, R) -integrals 



Representing the (^"-tensor by (9), one has the analogue of 
(16) for the (6, R) -integrals: 
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where the auxiliary vectors Q s read: 
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Since ( 0i ) 6 = 0, the analogue of the first term in (16) with 

vector Qg = Ei=i<2f (o°) 6 /(o) 6 vanishes - Equation (50) is 
well-known [9, 14, 16, 5]. We only mention it in order to show 
how nicely it fits into our general scheme. 

6. Discussion of the auxiliary vectors 

Finally we want to investigate properties of the auxiliary 
vectors Qq , defined in (7) . Of particular practical interest 
is the contraction of a tensor integral with a chord. 5 Usually, 
a scalar product qtk is expressed in terms of the difference of 
two propagators, which can then be cancelled such that a tensor 
integral of lower degree is obtained. This simplification is often 
the first step considering the original diagram, where one may 
further select q n = 0. Our approach offers an alternative due 
to the fact that the contraction of one of the vectors (7) with 
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a chord yields a simple expression (utilizing here the choice 

<?„ = 0): 
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Further we observe that the contraction of any of the extra terms 
introduced in (1 1) and (12) with a chord vanishes (for q n = 0). 
In this way, we have not only reduced the rank of the tensor by 
one, but we additionally obtain a simple expression in terms of 
the lower rank tensors. 

Relations (52) and (53) also allow a very simple derivation 
of (16) for q$ = by projecting with a complete set of chords. 
Writing (16) as: 
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with kqk = — 5 [cfc — c 5 + Yfe 5 — F55], from (52) and (53) im- 
mediately follows the equality. 

For 55 y£ a shift, k — > fe + (75, is needed on the integration 
momentum fc of the original integral. Under such a shift all 
vectors qi are shifted to ^ — > % — (75. The Q s for s = 1 . . . 5 
stay invariant, and the Qo shifts like a chord: 
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This is due to the fact that the Gram determinants are invariant 
under any shift and 
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After the shift one has only integrals with 175 = 0. Collecting 
all contributions and shifting back, one obtains (16). 



The auxiliary vectors Q s introduce in the recursions inverse 
powers of ()„. Relations (52) and (53) yield: 
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Equations (59) and (60) might suggest that after contractions 
with chords the effective inverse powers of ()„ are reduced. 
This is not true in general. For n = 5, however, this can be 
shown as follows: We write (52) as a system of linear equations 
for the vector Qq : 
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Xi = --{Y i5 -Y 55 ). 
Obviously, the matrix A satisfies 

A in g^A jv = q,qj, i, j = l,---4, 
and hence, 

det{q iqj ) = -i() B = -dot 2 (A). 
Therefore, 

A- 1 - () B -'- 
Solving the system (61) yields 
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— 1/2 

Thus the vector Qo is proportional to () 5 , while its additive 
components in (7) are proportional to O^ 1 . The same obviously 
applies for the auxiliary vectors Q s (s = 1 . . . 5). 

Summary 

We have presented a new, recursive reduction scheme for 
one-loop n-point tensor Feynman integrals, and derived explicit 
expressions covering tensors up to rank R = n, with n < 6. 
The crucial point is the derivation of relation (16), which ex- 
presses (5, R) -integrals in terms of (5,i? — 1)- and (4, R — 
1) -integrals, and of relations like (31), which express (k,R)- 
integrals (k < 4) in terms of (k, R — 1)- and (k — 1, R — 1)- 
integrals plus additional terms. Both types of representations 
make use of auxiliary vectors Q%, 

The recursive scheme is very convenient for explicit cal- 
culations and contains the complete calculational chain of ten- 
sor reduction, for both massive and massless propagators, and 
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works with dimensional regularization. The only necessary pack- 
age to be added is one for the evaluation of 1 -point to 4-point 
scalar integrals. For this we chose the library of scalar func- 
tions QCDloop/FF [26, 12] and made careful numerical checks 
of the recursive reductions. The numerical output was com- 
pared, whenever possible, to another, independent Fortran pro- 
gram [5, 10]. Additionally, 6-point functions have been com- 
pared with a program of P. Uwer for massive particles [27], and 
5- and 6-point functions with the Golem package [13], which 
applies for massless particles only. The highest rank covered 
was R — 5. We have agreement in all cases, both for the diver- 
gent and finite parts of the tensor integrals. Various tables with 
numerical comparisons may be found in [28]. 

Differing from the reductions derived in [5, 10], one can- 
not avoid the appearance of inverse Gram determinants ()s in 
the recursions. There is a lengthy discussion in the literature, 
and the rising numerical experience with explicit calculations 
of cross-sections leads to the conclusion that it is a small frac- 
tion of phase-space points which is substantially concerned. For 
those points, we have to switch to our code based on [5, 10], 
which is more appropriate to the specific kinematics around 
05 = 0. But doing so, the recursive character of the scheme is 
given up. It is also possible, in the recursion described here, to 
re-introduce and/or keep integrals in higher dimension in order 
to avoid inverse Gram determinants, as exemplified in another 
approach [14]. 
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